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The hadron inclusive spectra in pp collisions at high energies are analyzed within 
the soft QCD model, namely the quark-gluon string model. In addition to the sea 
quark distribution in the incoming proton we consider also the unintegrated gluon 
distribution that has an increasing behaviour when the gluon transverse momen- 
tum grows. It leads to an increase of the inclusive spectra of hadrons and their 
multiplicity in the central rapidity region of pp collision at LHC energies. 



1 Introduction 

As is well known, hard processes involving incoming protons, such as the deep-inelastic 
lepton-proton scattering (DIS), are described using the scale-dependent parton distri- 
bution functions (PDFs). Usually such distribution is parametrised as a function of 
the longitudinal momentum fraction x and the square of the four-momentum transfer 
q 2 = —Q 2 , integrated over the parton transverse momentum k t . However, for semi- 
inclusive processes, such as the inclusive jet production in DIS [I], electroweak boson 
production [T], etc., the parton distributions unintegrated over kt are more appropriate. 
The theoretical analysis of the unintegrated quark and gluon PDFs is presented recently 
in [2]. According to [2J, the gluon distribution function g(k t ) as a function of k t at fixed 
Q 2 has a very interesting behaviour at small x < 0.01, it increases very fast starting 
from almost zero values at kt ~ 0. In some sense, g(k t ) blows up when kt increases, 
then, it falls down at k t close to 100 Gev/c. In contrast to that the quark distribution 
q{kt), as a function of k t , is almost constant in the whole region of k t up to k t ~ 100 
GeV/c and smaller than g(k t ). These parametrisations of the PDFs were obtained in 
[2] within the leading order (LO) and the next to leading order (NLO) approximations 
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of QCD at Q 2 = 10 2 (GeV/c) 2 and Q 2 = 10 4 (GeV/c) 2 from known (DGLAP-evolved 
[3]) parton densities determined from global data analysis. At small values of Q 2 the 
nonperturbative effects should be considered for properly parametrising the PDFs. The 
nonperturbative effects can arise from the complex structure of the QCD vacuum. The 
instantons are one of the well studied topological fluctuations of the vacuum gluon fields, 
see, for example, [I]- [6] and references therein. In particular, it is shown [6] that the 
inclusion of the instantons results in the anomalous chromomagnetic quark-gluon in- 
teraction (ACQGI) which for the massive quarks gives the spin-flip part of it. Within 
this approach the very fast increase of the unintegratcd gluon distribution function at 
< kt < 0.5 Gev/c and Q 2 = 1 Gcv/c is also shown. 

2 Inclusive spectra of hadrons in pp collisions 

Let us analyze the hadron production in pp collisions within the quark-gluon string 
model (QGSM) [7] or the dual parton model (DPM) [8] including the transverse motion 
of quarks and diquarks in colliding protons [9l [10] . As is known, the cylinder type graphs 
presented in Fig.l dominate at high-energy hadronic interactions [7]. A physical meaning 
of the graph presented in Fig.l is the following. The left-hand side diagram of Fig.l, 
the so called one-cylinder graph, corresponds to the case when two colorless strings are 
formed between the quark/diquark {q/qq) and the diquark/quark (qq/q) of the colliding 
protons, then, at their breaking the quark- ant iquark and diquark-antidiquark pairs are 
created and fragmented to hadrons. The right hand-side diagram of Fig.l, the so called 
multi-cylinder graph, corresponds to a creation of the same two colorless strings stretched 
between valence quarks and diquarks and many strings stretched between sea quarks and 
antiquarks in the different protons. 




Figure 1: The one-cylinder graph (left) and the multi-cylinder graph (right) for the 
inclusive pp — > hX process. 
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The inclusive spectrum of hadrons in the QGSM is written as follows [51 [TP]: 

, da 
'tPp 



d °° 

p(x,p t ) = E-^r = ^2 o n (s)(j) n (x,pt) , (1) 



n=l 

where E, p are the energy and three-momentum of the produced hadron h in the c.m.s. 
of colliding protons, ^fs is the initial energy; x,pt are the Feynman variable and the 
transverse momentum of h respectively; a n is the cross section of n cut-Pomcron inter- 
action which is calculated within the "quasi-cikonal approximation" the function 
4> n (x,p t ) has the following form [5] : 

4>n(x,p t )= / dx\ / dx 2 tp n (x,p t ;xi,x 2 ) , (2) 

where 

iP n {x,p t ;x u x 2 ) = fW {x+ , p t ; Xl )F<£> (x_,p t ; x 2 )/F^ ) (0, p t ) + (3) 
+Fg> (x+ , Pt ; X! ) (x-,p t -,x 2 )/Fg))(p,p t ) + 
2(n-l)F(f{x + ,p t ;x 1 )F^\x.,p t ;x 2 )/F^){0, Pt ) . 



and x± = 0.5( \/x r +x\ ±x),x t = 2^J (m\ + pfj/s, 

^(aJi.pt^La) = I d 2 k t j^\x ±1 k t )Gr^h ( — ,h;p t ) 



(4) 



Here r stands for the flavour of quarks and diquarks, jv (x 1 ,kt) is the quark distribu- 
tion function that depends on the longitudinal momentum fraction x' and the transverse 
momentum k t . G T ^ h (z,k t ;p t ) = zD T ^ h (z,k t ;p t ), where D T ^. h (z,k t ;p t ) is the fragmen- 
tation function (FF) of a quark (antiquark) or diquark of flavour r into a hadron h. We 
assume that both the distribution and the fragmentation functions arc factorized over 
the longitudinal and transverse momentum. If their dependence on the transverse mo- 
mentum has a Gaussian form, then the function iv (x±,pt; x) can be written as follows 
0: 

F^(x ± ,p t ;x lt2 ) = f^(x h2 )G T ^ h (z)I n (z, Pt ) , (5) 
where the function I n (z, P t) reads [9] 

I n {z,p t ) = — exp(-7 z p^) 5 7z = __Z — p=~ . (6) 
7r 1 + npz z 7 
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Note that at x = the function Ft does not depend on n. In our calculations we have 

parametrized the quark distribution and FF by an exponent. In this case the function 

(n) 

I n has a more complicated form. However, in this case Ft does not depend on n too. 

• Gluon distribution in proton . 

As is mentioned in the Introduction, the unintcgrated gluon distribution in the proton 
at small values of function of kt, increases very fast when kt increases and then 

slowly decreases, according to the instanton vacuum approach for the massive quarks [6] 
at low values of Q 2 ~ 1 GcV/c. The similar behaviour for g(kt) was obtained within the 
NLO QCD calculations at large Q 2 = 10 2 GeV/c and Q 2 = 10 4 GcV/c for the massless 
quarks [2] . This allows us to assume that at Q 2 = and x = the gluon distribution in 
proton, as a function of kt, has a similar behaviour. 

g(k t ) ~ k^ eM-bgh) , (7) 

where the parameters a g > and b g > 0. 

• Hadron production in central rapidity region. 

According to the Abramovskiy-Gribov-Kancheli cutting rules (AGK) [12], at mid-rapiditiy 
only Mueller-Kancheli type diagrams contribute to the inclusive spectrum of hadrons. In 
our approach the function Ft is parametrized on such a way, that in the central region 
(y = 0), when x ~ and z ~ 0, it becomes proportional to n and satisfies to the AGK 
cancellation. Thus, 

oo 

Pq(x = 0,p t ) = <f>q(0,Pt) ^2 n(T n{s) = gS A (/) q (0, p t ) , (8) 
n=l 

where 4> q (x = 0,pt), depends only on p t . Considering the gluons from incoming protons, 
which may split into qq pairs, we get an additional contribution to the spectrum: 

oo 

p g (x = 0,p t ) = <pg(0,p t ) ^2(n- l)cr„(s) = (f> g (0,pt)(gs A - a nd ) , (9) 

n=2 

Where A = 0.12, g=21 mb and a nc i is the cross section of any number cut-Pomcron 
exchange. The quantities 

oo oo 

^ na n (s) = gs A ; ^ a n (s) = a nd (10) 

n— 1 n—1 

were calculated in |llj within the "quasi-cikonal" approximation [11| . Assuming that one 
of the cut-Pomerans is always scratched between valence quarks and diquarks which are 
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not coming from the splitting of gluons, in Eq. (9) we excluded unity from n. Finally, 
we can present the inclusive spectrum at x ~ in the following form: 

p(Pt) = Pq{x = Q,pt) + p g (x = 0,p t ) = gs A cf) q (0,pt) + (gs A - a n d) </> g (0,Pt) , (11) 

We fix these contributions from data on the charged particles pt distribution, parametris- 
ing them as follows: 

<t> q {0,Pt) = A q exp(~b q p t ) , (f> g {0,pt) = Ag^/jHexp(-b g pt). (12) 

The parameters are fixed from a fit to data on p t distribution of charged particles at 
y = 0: A q = 4.78 ± 0.16; b q = 7.24 ± 0.11 and A g = 1.42 ± 0.05; b g = 3.46 ± 0.02. 

3 Results and discussion 

In Fig. 2 we illustrate how our approach works for the description of the ISR data on 
inclusive spectra of charged pions and kaons. The solid lines correspond to our calcula- 
tions without the gluon contribution p g . As one can see, we have obtained a quite good 
description of the spectra as functions of p t , up to p t = 1A GcV/c at different values of 
x. Though, in Fig. 2 we present results for y/s = 53 GeV, but the description of data at 
other ISR energies (y/s = 23.3, 30.6, 44.6, 53 GeV) is similarly good. The result of the fit 
to data on the charged hadron inclusive spectra is presented in Fig. 3. The long dashed 
curve corresponds to the quark contribution p q (x = 0,pi) given by Eq.(j8|), whereas the 
short dashed line is the gluon contribution p g (x = 0,pt) (Eq.©) to the invariant yield 
d 3 N/dyd 2 pt; the solid curve corresponds to the sum of both contributions, see Eq. (|T2|) . 
One can see that the conventional quark contribution $ l? (j/ = 0,Pt), see Eq. (|T2|) . is able 
to describe the data up to <1 GcV/c, whereas the inclusion of the gluon contribution 
allows us to extend the range of good description up to 2 GeV/c. At larger values of pt 
the contribution of hard processes is not negligible and one has to take them into account 
based on the perturbative QCD. 

4 Conclusion 

Our study has shown that the soft QCD or the QGSM is able to describe the experi- 
mental data on inclusive spectra of light hadrons like pions and kaons produced in pp 
collisions at not large values of the transverse momenta p t < lGeV/c rather satisfactorily. 
The inclusion of the unintegrated distributions of gluons in the proton due to instanton 
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vacuum excited in the strong pp interaction allows us to extend the satisfactory descrip- 
tion of the experimental data on the inclusive spectra of charged hadrons at y = up to 
p t ~ 2GeV/c. At higher values of the transverse momentum the hard parton interactions 
should be considered for describing the experimental data. 




Figure 2: The inclusive spectra of tt (left) and K (right) mesons in pp collision 
Eda/d 3 p [mb GcV~ 2 c 3 ] at ^fs =53 GeV compared with the ISR data [13]. The solid 
lines correspond to our calculations whithout the gluon contribution p g 
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Figure 3: The charged particles yield in the central rapidity region (y = 0) as a function 
oipt at -^=540 and 900 GeV (top) and y/s=1.8, 2.36 and 7 TeV (bottom) compared with 
UA1, CDF and CMS data [HJ [T5J [TB] . The long dashed curve corresponds to the quark 
contribution p q (x = 0,Pt) given by Eq.([8]), whereas the short dashed line is the gluon 
contribution p g (x = 0,pt) (Eq.©) to the inclusive yield d 3 N/dyd 2 p t ; the solid curve 
corresponds to the total calculation including both these contributions, see Eg . (fT2"l) . 



